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I. 


INTRODUCTION 


This technical memorandum presents a program for an HP-25 or HP-25C hand 
calculator that will calculate accurate LORAN-C time differences. The program is 
most useful when checking the accuracy of a LORAN-C receiver at a known latitude 
and longitude without the aid of an expensive computer. It can thus be used to compute 
time differences for known landmarks or waypoints to predict in advance the approximate 
readings during a navigation mission. 

The Fischer spheroid is the earth model used in the program. This is the 
same system the U,S, Coast Guard and the U, S, Defense Mapping Agency use to 
calculate transmitter baseline lengths and LORAN-C navigation charts and tables. 

It is assumed that the program operator is familiar with the general operation 
and programming capabilities of the HP-25 calculator 

II, PROGRAM DESCRIPTION AND ACCURACY 

The forty -three step program has a running time of approximately ten seconds 
and is relatively simple to operate. It operates in the decimal degree mode. Latitudes 
and longitudes must be stored as decimal degrees (dec. deg,). For each receiver position 
the program is run once to solve for T^/ Figure 1, and again to solve for Tj], Formula 
I, Figure 1, will solve for T]/ the time delay for the first master-slave pair. Run the 
program again to solve for slave Tj 2 again use formula I to solve for T 2 , the time 
delay for the second master-slave pair, etc. Storage register seven (7) is left open 
to hold whatever number or time delay is >ired. 

The indicated formulas in Figure 1 are omitted in this program due to a lock 
of program memory. However, this has only a minor effect on accuracy. The comparison 
with the computer values. Table 1, shows that the hand calculator is usually 
within ±5ps. Although the program has not been tested over a large geographic area, 
it has always been within ± 10 |Js of computer results compared so far. 


III. COMPUTING TIME DIFFERENCES 

By following the HP-25 program forms. Appendix A, time differences can be 
easily calculated. An example problem is shown on page 3 using Athens, Ohio coordinates 
in the supporting data. Table 1, and also using the U,S. East Coast Chain with the rate 
of 9930(SS7). The instruction steps in the example correspond to the instruction steps 
in Appendix A, Data for other LORAN chains is found in Appendix B 
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IV. FUTURE PROGRAM CONSIDERATIONS 


Similar programs for other makes of hand calculators, such as the TI-58, would 
be useful. Card reading calculators should hold the entire sequence of formulas, in- 
cluding those left out of this program. Card reading calculators should also be able to 
change time differences to latitude and longitude coordinates by using one of the 
methods described in Navigation It is possible that the HP-25 can perform one of 
these simple methods, but the accuracy and utility has yet to be determined. 


V. SUMMARY 

Although card reading calculators should give best results, the HP-25 LORAN 
Time Difference Program has proven itself to be reliable over the coordinates that it 
has been tested. Changing time differences into position on a low-cost hand calculator 
is a future area of pursuit. 
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HP-25 Program Form 


TitIr Exampl e Probl em Page. 

Programmer LOHAN Time Differe nces for Athens a t r^te 9930 (SS?) . 




INSTRUCTIONS 


Key in program _ 



Latitude, 
Longitude, 


Compute T 


Store T 


Latitude, <j) 


Longitude, Ar 


Return to top of memo 
Compute T 


Compute T - T 

a. JL ... 

Compute T 


INPUT 

DATA/UNITS 


21^32.339 

0.99664767 

39.1930 

82.0615 


34.034604 


77 . 3^76 



27.015849 


80.065352 




12 Compute another 


secondary, go to step 


Latitude, (pr 


Longitude, A-r 


Return to top of memor 


4602.0 


11000 


2695.57 


46.463218 


53.102816 



12 Compute another 

secondary, go to step 7 
chain, go to step 4; 
receiver position, go 
to step 3 
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The LORAN reading or Hme difference (T) is defined by the formula below and 
portrayed by Figure 1 . 

* Indicated formulas are omitted. A =21282.33935 psec. (Equotoriol rodius) 

(Minor geodetic corrections) B = 21210.99393 psec. (Polar radius) 

Formula (I) T = P + A + T - T C = B/A = 0.99664767 

s rn 


Master Slave 

Station p Station 



- ton ^ I C fan 1 

X = cos Msin(p|^)sin(p^)+cos(P|^)cos(P^)cos(X|^-X^)l 

* P^ = (A-B)l X-sin(X)l/(4.0)ll.0+cos(X)l 

* Ao=lsin(P|^) + sin(P^)l^ 

* Q = (A-B) 1 X +sin(X)l /(4.0) 1 1 .0-cos(X) 1 

* B^ =lsin(p|^)-sin(P^)] ^ 

d = AX-A P-B Q^AX 
o o 

d = T if (4>_,X )=(<p ,X ), Master 
m I I m m 

D = T if (0_, XJ =(d) , X ), Slave 
s II s s 

T = p +A+ T - T 
s m 


Figure 1. Method for Computing LORAN Hyperbolic Lattice. 
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Place 

Latitude 

Longitude 

Rate 


M-X 

M-Y 

M-Z 

Athens, OH. 
39® 19* 30" 

9930 

15971.836 

15979.0 

42354.859 

42353.1 

53738,113 

53734.8 

67801,313 

67797.9 

82® 06* 15" 

9960 

16591.750 

16593.2 

28539.574 

28538.3 

42594.543 

42596.5 

56775.641 
56774. 1 

Lamar, PA, 
41® 07* 00" 

9930 

16355.363 

16362.0 

40578.895 

40574.4 

52022.266 

52016.7 

68816.563 

68811.3 

77® 28* 00" 

9960 

16488.363 

16489.6 

28448.500 

28447.4 

44219.320 

44223.8 

59415.660 

59414.8 

Corning, N.Y, 
42® 10* 00" 

9930 

16361.043 

16367.4 

39928.367 

39923.0 

51536.559 

51530.4 

68586.563 

68580.9 

77® 06* 00" 

9960 

16559.109 

16560.5 

28759.621 

28759.4 

45016.137 

45022.0 

59982.484 

59982.6 

Logan, OH. 
39® 31' 30" 

9930 

15952.5 

15958.2 

42286.7 

42281.2 

53690.5 

53686.4 

67605.6 

67602.6 

82® 22* 55" 

9960 

16594.6 

16594.4 

28585.9 

28583.3 

42688.4 

42690.5 

56673.8 

56671.9 


Nof’el The compufer results are obtained by using the complete algorithm of 
Figure 1. 


Table 1, Supporting Data, HP-25 program results are listed below 
each computer result. 
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A. Program end ^rogram Form 


TiIIl*. L<^^1AN \ \N\t l>iri I KrM<'-L. P.i(je_/__of 2 

Svvildi to PRGM mot!.', pr.-.b » ! ' pmcm j , tlion kry i” the prfxjMrn. 


tiiriPLAY 

\.v r 


Y ' 

1 

7 

T 

CCdMMfcNTS 

HruiSIERS 

lINt 1 

coot 1 

‘NIHY 1 

^ i 




Oij j 
oi I 





i 

1 



Rq - 

M <\ 
U / 1 1 
Hj 

Hcl i‘i 

',•7 L 
( <’ 't >\ j 

X 

( 

! 

1 

1 

1 




1 

0^ 

'! t 

(. 

1 


^'c Pt_ 


0'^ 

C.oV'l r 
. ' '1 n<l t 

C 



R , lt>T 

1 

01 1 




(,<Ut.jl*.*u) 

1 » 





— 



IS C'l 

‘J tw\ 

f<C L 
RCL 7 

f tan 

X 

(j'r 

c 



1 



Uo 

f^’r 




R 2 A t 

07 

vLa 

t: 





oa 

jiCi 

i 1 

t en _ 
,tan ^ R_ 

C; 





09 

fr 





10 


a tan' 

fir 



Cem outing X 

((Lv-cLi|)L 

4 

11 


t C.C-S 

>pw 1 

■*> 

Dr 



.iiA 

12 

_2l 

X ‘ Y 1 

(5r 1 

C t 




R 4 A Ji. 

13 

11 7^ 

Las^ X. 

tf/p ! 


0’C^2k_ 



(.Lu*i1j£.3) 

14 

a 

{ Sin 

j 

- 

pit 

c^ji 



15 

2i 

X - y 





AVmjaiUYii^ ih-K.- 

R 5 .5.1212.- 

16 

lid 

•f s in 

binpT 


CC-^» - 


stack __ jQ\clixtaA*\!r. _ 

M'lll _ 

17 

li 71 

Last X 

fir 


Kh.l^x. 

cii5A 

an cEtn .rcAisltiLCT). 


18 

liCT 

•? CC5 

i OS/ft 

S/O/ft 






19 

22 

R ^ 

kir^ 

S 1 n . _ 

OOSijL- - 

C0'^5r 



20 

^ 1 

X 

S 

LCSfit 











S 


R, 


HKfl 

X 




s 




KM 

X ? > 

c<y^ 


^MWHt 

r 

J 



24 

22 

R *■ 

C6VvC'<34ff 


s 

C<-\3&: _ 



25 

Zi a 

ECL ^ 

X It 


S' 

s_ 



26 

2VC2 

Rcl 2 

A T 

X 8 

CCM^Coy;, 

S). - - 



27 

11 


iL-Xt 


S 

5_ 



28 

liCT 

£ccs 

Crs^XA'XT) 

C6SrxCC-VY 

5_ .. - 

S 



29 


X 


S 

S 

5.. __ 

|^= Co;:,S^Cc^LC^i[)^ 


30 

5i 

+• 

R+S 

s 

•S 

s 



31 

ISOS 

fl CCS * 


.s 

.$ 

s 




32 

Cl 

1 

1 

X' 

S 

s 

X. .f « 


33 

08 

6 

|3 

A - 

S 

. . . 

jlcarcki- t« r>uii\ub- 


34 

oc 

0 

l?C 

X . - 

S 

s 



35 

\S7? 

3 11 

iT* 

1?C 

s 

s 



36 

1 ( 


i?c/rr 

X 

■S. .. . 

s 



37 

7( 

4- 

X ('“‘lO 

s 

s 

s 



38 

2*105 

RCl 5 

A_. 

X C'faAl 

s 

•S 



39 

0 J 

X 


IS 

V 

s 

l4 lAastcr 0 


40 

7*1 

R/<f 

1 S 

s 

5 

IV Slaec, K/S _ 


41 

27 OC 

J^C L 0 

Tm 

1 5 





42 

U 


■TI_-T,n 






43 

13 OC 

&T0 cx 

Ts “ Tm 




AdJLA f ^ 9 


44 







XIrvl. X 


45 









46 









47 

" ' “ 








48 





j 




49 





1 


J 1 


At’ ^ 

of 


HEWl. FIT ^PACKARD 
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HP-25 Program Form 


Title LQRAN IimcLPjff®rence 

Programmer 


INSTRUCTIONS 


Key in program 
Store constants 


. Page_^__of 


Store the receiver 


coordinates. . _ 

LatiJude,_‘^R._ 

Longitude, 

Store the primary 

I'-iTismitter coordinate: 

Lot'» lo'e, <t>T 


Longitude, Xj 

Comp Jte .Ln 

Store T m 

Stwe the secon dary 

transmitter coordinate 
Latitude, Oj 


Longitude, 


Return to top of memory 


Compute 

__Com^te Tj-Tjjj 

Compute T by 

adding the coding 
delay (A) Baseline 
length p for the seconc 
To c^pute T for onother 

secondary, go to s tep 

chain, go to step 4; 

receiver ^osjtion, to 

to step 3r etc . 
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B. LORAN Chain Data 



U.S. Eait Coast Chain - Ratt 9930 (SS7) 


Station 

Coordinates 

Latitude & 

Longitude' 

Station 

Function 

Coding 
Delay & 
Baseline 
Length 

Radiated 

Peak 

Power 

Carolina 
Ileach. NC 

34 03-46.04 N 
77-5446.76W 

Master 

... 

700 kW 

J'lpitcr, 

Florida 

2701-58.40N 

80-06-53.52W 

VI 

Secondary 

11.000;is 
2695. 57ps 

300 kW 

Cape Race. 
Newfoundland 

46 '46-32. 18N 
53-1 0-28. 16W 

X 

Secondary 

28.000ms 

8389.55ms 

1.8MW 

Nantucket. 

Massachusetts 

41-15-11.93N 

69-53-39.09W 

Y 

Secondary 

49.000ms 

3541.27ms 

300 kW 

Dana. 

Indiana 

39 51-07.54N 
87-29-12.14W 

Z 

Secondary 

65.000ms 

3560.68M< 

400 kW 


For Jupiter, Florida 
Coding delay A = 1 1000 ps 
Baseline P =2695,57 ps 


I 


Northeast U.S.A. Chain (reconfigured 7/78) » Rats 9960 (SS4) 


Station 

Approximate 

Latitude & 
Longitude 

Station 

Function 

Coding 
Delay ft 
Baseline 
Lertgth 

Radiated 

Peak 

Power 

Seneca. 
New York 

42-42-30. '6C 
76-49-33 

Master 

... 

1.0 MW 

Caribou. 

Maine 

4648-27.20N 

67-55-37.71W 

W 

Secondary 

11,000ms 

2797.20 

350 KW 

Na.i bucket. 
Massachusetts 

4M5-11.93N 

6958-39.09W 

X 

Secondary 

25.000ms 

1969.93 

300 KW 

Carolina Beach 
N. Carolina 

34-03-46.04N 

77-54-46.76W 

Y 

Secondary 

39.000MS 

3221.63 

700 KW 

Dana. Indiana 

3951 0/.54N 
87-29 12.14V7 

Z 

Secondary 

34000 

3162.06 

400 KW 
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